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ABSTRACT 

Using mean-field models with a dynamical quenching formalism we show that in finite do- 
mains magnetic helicity fluxes associated with small-scale magnetic fields are able to alleviate 
catastrophic quenching. We consider fluxes that result either from advection by a mean flow, 
the turbulent mixing down the gradient of mean small-scale magnetic helicity concentration, 
or the explicit removal which may be associated with the effects of coronal mass ejections 
in the Sun. In the absence of shear, all the small-scale magnetic helicity fluxes are found to 
be equally strong both for large-scale and small-scale fields. In the presence of shear there is 
also an additional magnetic helicity flux associated with the mean field, but this flux does not 
alleviate catastrophic quenching. Outside the dynamo-active region there are neither sources 
nor sinks of magnetic helicity, so in a steady state this flux must be constant. It is shown 
that unphysical behavior emerges if the small-scale magnetic helicity flux is forced to vanish 
within the computational domain. 
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1 INTRODUCTION 

Both mean-field theories as well as direct simulations of the gen- 
eration of large-scale magnetic fields in astrophysical bodies such 
as the Sun or the Galaxy invoke the effects of twist. Twist is typ- 
ically the result of the Coriolis force acting on ascending or de- 
scending magnetic field structures in a stratified medium. The net 
effect of this systematic twisting motion on the magnetic field 
is called the a effect. In text books the a effect is normally in- 
troduced as a result of helical turbulence ( lMoffatt|[T973 : |Parkeil 
Il979l : iKrause & Radler|[T980l), but it could also arise from mag- 



netic buoyancy instabilities (*Sc hmittll987l . [Brandenburg & Schmitj 
1998) . The latter may also be at the heart of what is known 
as the Babcock-Leighton mechanism that describes the net ef- 
fect of the tilt of decaying active regions. Mathem atically, thi s 
mechanism can also be described by an a effect ISti}3 Il974h . 
Regardless of all these details, any of these processes face a 
serious c hallenge connected with the conservation of magnetic 
helici t v jPouguet. Frisch & LeoratI 1 19761: iKleeorin & Ruzmaikii3 



1 19821 : iKleeorin. Rogachevskii & Ruzmaikirj Il995t) . The serious- 
ness of this is not generjilly appreciated, even though the conser- 
vation of magnetic helicity has lon g been associated with what 
is called ca tastrophic a quenching jOruzinov & Diamond 1 19941 
1 199511 19961) . Catastrophic a quenching refers to the fact that the 
Q effect in helical turbulence in a periodic box decreases with 
increasing magnet ic Reynolds num ber for equi partition stren gth 
magn etic fields dVainshtein & Cattanedl 19921 : ICattaneo & HughesI 
1 19961) . This would be 'catastrophic' because the magnetic Reynolds 
number is large (10^ in the Sun and 10^'' in the Galaxy). 

A promising theory for modeling catastrophic a quenching 



in a mean-field simulation is the dynamical quenching formula, 
i.e. an evolution equation fo r the a effect t hat follows from mag- 
netic helicity conservatio n l lKleeorin & Ru zmaikin 1982). Later, 
iField & Blackmad ( l2002h showed for the first time that this for- 
malism is also able to describe th e slow saturation of a helical dy- 
namo in a triply-periodic domain ( lBrandenbur32001iJ) . As this dy- 
namo runs into saturation, a large-scale magnetic field builds up, 
but this field possesses magnetic helicity. Indeed, the eigenfunc- 
tion of a homogeneous dynamo has magnetic and current he- 
licities proportional to a. However, this concerns only the mean 
field, and since the helicity of the total field is conserved, the small- 
scale or f luctuating fie ld must have magnetic helicity of the oppo- 
site sign I Seehafeill 1996.). This l eads to a reduction of the a effect 
dPouQuet. Frisch & LeoraJl 19761) . 

The dynamical quenching formalism is now frequently used 
to model the n onlinear behavior of mean- field dynamos with and 
without she ar jBlackman & Brandenbur i|2002), open or closed 
boundaries terandenburg & Subramanian' 2005), an d sometimes 
even without a effect I Yousef, Brandenburg & Rti digeil 1200 3l : 
Brandenbu rg & Subramaniani [2005). However, it became soon 
clear that the catastrophic quenching of the a effect can only 
be alleviated i n the presence of magne ti c helicity fluxes out 
of the domain jBlackman & Field l2000allbl IKleeorin et alj|200d. 
[2002). T here are various coritrib u tioris to the magnetic helic- 
ity fl ux teogachevski i & Kleeorin 2000; ' Vishniac & Choi 1200 ll : 
.Subramanian & Brandenburg.2004..2006. ), but one of the most ob- 
vious ones is that associated with advection. lShukurov et al. U200^ 
have implemented this effect in a mean-field model with dynamical 
quenching in order to model the effects of a wind on the evolution 
of the galactic magnetic field. One goal of the present paper is to 
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study this effect in more detail. In particular, it is important to clar- 
ify the consequences of boundary conditions on the local dynamics 
away from the boundaries. Indeed, is it really true that a helicity 
flux has to be maintained all the way to the boundaries, or can the 
helicity flux be confined to part of the domain to alleviate catas- 
trophic a quenching at least locally? What happens if this is not 
the case? 

The notion of alleviating catastrophic a quenching only 
locally is sometimes invoked in models of the solar dynamo 
that rely on the production of strong magnetic fields at the 
bottom of the convection zone. By placing the a effe ct only 
near t he surface, as is done in the interface dynamo of IParkeJ 
( Il993[) or dynamos that are controlled b y meridional circula- 
tion jChoudhuri. Schiissler & Dikpatilll995 l). one may evade catas- 
trophic quenching more easily. On th e other hand, as shown by 
lYousef. Brandenburg & RudigeiH200 3'), the effects of magnetic he- 
licity conservation can even play a role even if there is originally no 
a effect. It is therefore important to understand in more detail the 
physics of dynamical a quenching and its dependence on magnetic 

helicity fluxes. 

Our starting point in this paper is the model of lShukurov et al.l 

( I2006") . where magnetic helicity fluxes were driven by the advection 
from a wind. This allows us to study the effects of varying strength 
of this flux in different parts of the domain. For simplicity, and in 
order to isolate the main effects, we ignore shear in most parts of 
this paper. In view of later applications to the Sun and the Galaxy 
this is clearly artificial, but it helps significantly in the interpretation 
of the results. In particular, in the absence of shear, it is possible to 
have steady solutions, or at least solutions whose magnetic energy 
density is constant in time. This simplifies the interpretation of the 
results. 



E = V{SxAy) - SAyX + U X B + £ - 77^0 J, 



(2) 



where U is now the flow associated with the wind only and does 
not include the shear flow. Next, we express B = V x A in 
terms of the magnetic vector potential A, and solve equation ^ 
in its uncurled form, dA/dt = —E — 'V(j), where (j) is the 
mean electrostatic potential. We perform a gauge transformation, 
A ^ A + VA, with the choice A = J (0 ~ 
removes the gradient term to yield 



SxAy) dt, which 



OA 
'dt 



(3) 



which is then the final form of our equation for A. This form of the 
equation together with boundary conditions for A characterize the 
gauge used to calculate magnetic helicity densities and magnetic 
helicity fluxes for the mean field. 

We solve equation l[3j in the domain < z < L and as- 
sume either a vacuum or a perfect conductor boundary condition 
on 2: = L. This means that on z = L the mean magnetic field ei- 
ther vanishes, i.e. ~ By = 0, or that its z derivative vanishes, 
i.e. Bx.z ~ By^z — 0, where a comma denotes partial differentia- 
tion. In terms of A this means that on 2; = L we have either 



Ax,Z Ay 



(vacuum condition), 



Ax = Ay = (perfect conductor condition). 



(4) 



(5) 



It is well known that the solutions can be in one of two pure 
parity states that are either symmetric (S) o r antisymmetric (A) 
about the midplane l lKrause &Radleil[T980h . so we have either 
Bx,z = By^z = or Bx = By 
this means either 



on 2; = 0. In terms of A 



2 THE MODEL 

2.1 Evolution equation of the mean-field 

In this paper we consider a simple mean-field dynamo in a local 
one-dimensional domain. Such a model could be applicable to one 
hemisphere of a rotating disc or to the region close to the equa- 
tor of outer stellar convection zones. Denoting the mean magnetic 
field by B = B{z,t), the coordinate z would correspond either 
to the height above the midplane in the case of the disc, or to the 
latitudinal distance from the equator in the case of a spherical shell. 
The X and y components would correspond to poloidal and toroidal 
fields, although in the absence of shear the two are interchangeable 
and cannot be distinguished. Using V ■ B = dBz/dz = 0, we 
have Bz ~ const = 0, i.e. no Bz field is imposed. Such a mean 
field could be obtained by averaging the actual magnetic field over 
the X and y directions on a Cartesian domain. 

The evolution of B is governed by the Faraday equation 



Ax = Ay = on z = (S solution) 



OB 

dt 



= -V X 



(1) 



where iS = —(C7s + x B — f + 77/10./ is the mean electric field, 
U is the mean flow in the z direction. Us ~ (0, Sx, 0) is a linear 
shear flow, £ is the mean electromotive force, J = V x B/jio 
is the mean current density, and /iq is the vacuum permeability. In 
one case we adopt a shear parameter S that is different from zero. 
Since the shear is linear, we can write Us x B as —SAyX plus a 
gradient term that can be removed by a gauge transformation. Thus, 
we have 



Ax,: 



Ay 



on z = (A solution). 



(6) 



(7) 



We note that the particular boundary conditions (|5} and ^ fix the 
value of A on z = 1/ or 2: = 0, respectively. In all other combi- 
nations the value of A is no t fixed and the magnetic helicity could 
exhibit an unphysical drift dBrandenburg. Dobler. & SubramanianI 
2002). However, in the present paper we study magnetic helicity 
density and its flux only in situations where either ([5j or l[6ll are 
used. 

We recall that, even though there is no effect, i.e. no mean 
flow in the y direction, we shall allow for a flow U in the z direc- 
tion. In a disc this would correspond to a vertical wind, while in a 
star this might locally be associated with meridional circulation. 



2.2 Magnetic helicity conservation 

In this paper we will study the evolution of magnetic helicity of 
mean and fluctuating fields. In our gauge, the evolution of the mag- 
netic helicity density of the mean field, Hm — A ■ B, is given 
by 



dh 



— = 2£-B- 2-qiioJ B V Fn„ 
ot 



(8) 



where i^m = E x A is the flux of magnetic helicity of the 
mean magnetic field. Under the assumption of scale separation. 
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ISubramanian & Bra ndenbur ? ('2006') have defined a magnetic he- 
licity density of the small-scale field in terms of its mutual link- 
ages. They derived an evolution equation for the magnetic helicity 
density of the small-scale field, 



dht 
dt 



-2£ ■ B - 2ri^ioj - b-V ■ Ft, 



(9) 



which is similar to equation l[8j, except that the £ ■ B appears with 
the opposite sign. This implies that turbulent amplification and dif- 
fusion of mean magnetic field (characterized by the £ term) cannot 
change the total magnetic helicity density, h = /im + hi, which 
therefore obeys the equation 



^ = -2r,^o J B V 
at 



(10) 



where F — F^i + F[ is the total magnetic helicity flux, and 
JB — JB+jb is the total current helicity density. 



2.3 Dynamical quenching formalism 

In order to satisfy the evolution equation for the total magnetic he- 
licity density l llOt . we have to solve equation (|9} along with equa- 
tion l[3), which implies that equations ([8} and dlOb are automatically 
obeyed. We make the assumption that the turbulence is at small 
scales nearly isotropic. This means that /loj ■ b ~ kfh{. The j ■ b 
term a lso modifies the mean electromo tive force by producing an a 
effect dPouquet. Frisch & Leoratlll976h . This is sometimes referred 
to as the magnetic a effect. 



^-rj ■ b/p. 



(11) 



where r is the correlation time of the turbulence. In the following 
we ignore compressibility effects and assume that the mean density 
p is constant Next, we assume that the turbulence is helical, so 
there is also a kinetic a effect proportional to the kinetic helicity. 



ax 



(12) 



where u; = V x u is the vorticity. The total a effect is then 

a — UK + au, (13) 

and the resulting mean electromotive force is 



£ = aB — ijtfioJ, 
where 



(14) 



(15) 



is the turbulent magnetic diffusivity. In the following we consider 
Tjt and r] as given and define their ratio as the magnetic Reynolds 
number, Rm ~ Vt/v- We shall express the strength of the magnetic 
field in terms of the equipartition value, 



(16) 



which allows us to determine r in the mean-field model via ir = 
PoT^Vt I Bcq ■ We characterize the value of rj in terms of the magnetic 
Reynolds number, Rm — rji/rj. With these preparations we can 
write the dynamical quenching formula as 



dt 



-2ritk[ 



£ B 



+ 



Rm 



OZ 



(17) 



^ Note that a constant mean density implies that there must exist a small- 
scale mass flux compensating the losses associated with the mass flux 'pU. 



where is related to the mean magnetic helicity flux of the fluc- 
tuating field via 



^ a 



poprjtkf — 



-'-'cq 



Fc. 



(18) 



In order to compute mean-field models we have to solve equa- 
tion l[3j together with equation l ll7b using a closed expression for 
the flux J- a- In this paper we focus on the advective flux propor- 
tional to auU , but in some cases we consider instead the effects 
of a turbulent magnetic helicity flux that we model by a Fickian 
diffusion term proportional to — Vqm, where Hc is a diffusion 
term that is either zero or otherwise a small fraction of r]t . 

In addition, we consider cases where we model magnetic 
helicity fluxes by an explicit removal of hi from the domain 
in regular time intervals At. Such an explicit removal of mag- 
netic helicity associated with the fluctuating field may model 
the effects of coronal mass ejections, although one would ex- 
pect that in reality such an approach also implies some loss of 
magnetic helicity associated with the large-scale field. The ap- 
proach of removing the fluctuating magnetic fi eld was employed by 
ISrandenburg. Dobler. & Subramaniag ( |2002| ) in connection with 
three-dimensional turbulence simulations to demonstrate that it is, 
at least in principle, possible to alleviate catastrophic quenching by 
an artificial filtering out of small-scale turbulent magnetic fields. In 
the present paper we model the occasional removal of hi by reset- 
ting its values 



hi 



hi — A /if in regular intervals At, 



(19) 



where Ahi = ehi is chosen to be a certain fraction e of the cur- 
rent value of hi. In our one-dimensional model the corresponding 
expression for the flux AFf can be obtained by integration, i.e. 



AFi{z,t) = / Ahi{z',t)dz'. 



(20) 



Since magnetic helicity densities and their fluxes are proportional 
to each other, we have simply 

J- a = OImU — Ha -7^ h AT" a , 

OZ 



(21) 



where AJ^a ~ {po7>V-t^f / Bcq)AFi is defined analogously to 
equation dlSl l. 

We note that the a effect will produce magnetic fields that 
have magnetic helicity with the same sign as that of a, and the 
rate of magnetic helicity production is proportional to aB^. In the 
northern hemisphere we have a > 0, so the mean field should have 
positive magnetic helicity. We recall that shear does not contribute 
to magnetic helicity production, because the negative electric field 
associated with the shear flow. Us x B, gives no contribution 
to magnetic helicity production, which is proportional to E ■ B, 
but it can still give a contribution to the flux of magnetic helicity. 
This is also evident if we write shear using the —SAyX term in 
equation after multiplying with B and using = dAy/dz, 
we find that this term can be integrated to give just an additional 
flux term, | SAy . However, this contribution belongs clearly to the 
magnetic helicity flux associated with the large-scale field and is 
therefore unable to alleviate catastrophic quenching. 



2.4 Model profiles and boundary conditions 

We consider a model similar to that of lShukurov et al.l ( l2006h who 
adopted linear profiles for aK and U of the form ax = aoz/H 
and U z — Uqz/H, where the height H was chosen to be equal to 
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Figure 1. Profiles of a and U for Waki = 0.2 and wuki = 1. 



the domain size, H — L. However, in order to separate boundary 
effects from effects of the dynamo we also consider the case where 
we extend the domain in the z direction and choose L — AH and let 
uk go smoothly to zero at 2 = and U ^ either goes to a constant 
for z > H or 'ii also goes smoothly to zero. Thus, we choose 



a = ao — Q{z;H,Wa), 
M 

where we have defined the profile function 

z-H' 



e{z;H,w) 



1 — tanh ■ 



(22) 



(23) 



which is unity for z <ti H and zero otherwise, and w quantifies the 
width of this transition. For the wind we choose the function 



u. = u,jj[i + {z/Hr 



-1/n 



9(z; Hu, wu), 



(24) 



with n = 20. Both profile s are s hown in Fig.[T] The strictly linear 
profiles o f IShukurov et^ ( l2006h can be recovered by taking L = 
H, Wa —> 0, and n —> 00. 

As length unit we take ki = n/2H, and a s time unit we take 



This deviates from IShukurov et al. who used 



n/H as their basic wavenumber. Our motivation for this change 
is that now the turbulent decay rate is equal to rjtki, without an ex- 
tra 1/4 factor. We adopt nondimensional measures for ao, Uo, and 
S, by defining 



(25) 



Cu = and Cs = —rj- 

To match the parameters of IShukurov et al.l i2006t) . we note that 
Cu = 0.6 corresponds to their value of 0.3, and the value kf /k\ = 
10 corresponds to their value of 5. 

We obtain solutions numerically using two different codes. 
One code uses an explicit third-order Runge-Kutta time stepping 
scheme and the other one a semi-implicit scheme. Both schemes 
employ a second order finite d ifferences. We begin b y reporting re- 
sults for the original profile of lShukurov et al.l ( l2006i) with L = H. 
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Figure 2. Space-time diagrams for and By for the marginal values of 
Ca for L = H with Cjj = and either the symmetric solution (S) with 
a vacuum boundary condition on z = H or the antisymmetric solution (A) 
with the perfect conductor boundaiy condition. In both cases the critical 
value Ca = 5.13 is applied. Light (yellow) shades indicate positive values 
and dark (blue) shades indicate negative values. 



3 RESULTS 

3.1 Kinematic behavior of the solutions 

When the magnetic field is weak, the backreaction via the Lorentz 
force and hence the om term are negligible. The value of Rm 
does then not enter into the theory. The effects of magnetic he- 
licity fluxes are therefore not important, so we begin by neglecting 
the wind or other transporters of magnetic helicity. For the linear a 
profile we find that the critical value of Ca for dynamo action to 
occur is about 5.13. These solutions are oscillatory with a dimen- 
sionless frequency uj = uj/rjtkf — 1.64. The oscillations are asso- 
ciated with a migration in the positive z direction. This is shown in 
Fig. |2] where we compare with the case of a perfectly conducting 
boundary condition at z = H for which we find C"'' — 7.12 and 
= 2.28. 

The fact that there are oscillatory solutions to the dynamo 
is perhaps somewhat unusual, but it is here related to the fact 
that a changes sign about the equator. Similar behavior has been 
seen in some other dynamos where a changes sign with depth 
(Stefani & G erbeth 20 03; R udiger. Elstner & Ossendrjjver 20()J; 
iRudiger & HollerbachI |2004| ; loiesecke', Ziegler & Rudigerl fe'oO^ 
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and in simulations of helically forced turbulence with a change of 
sign about the equator ( Mitra et al. 2009). In the latter case, how- 
ever, the outer boundaries were perfectly conducting. In our mean- 
field model such a case is also oscillatory, as will be discussed be- 
low. 

Note that we have made here the assumption that the solutions 
are symmetric about the midplane, i.e. Bi{z,t) — Bi{~z,t) for 
i = X or y. For the application to real systems such a symmetry 
condition can only be justified if the symmetric solution is more 
easily excited than the antisymmetric one for which Bi{z,t) = 
—Bi{—z, t) fox i = X or y. This is indeed the case when we adopt 
the vacuum condition at z = H, because the antisymmetric solu- 
tion has C"'* = 7.14 in that case. However, this is not the case 
for the perfect conductor boundary condition for which the anti- 



symmetric solution has Ca 



5.12. We remark that there is a 



striking correspondence in the critical Ca values between the anti- 
symmetric solution with perfect conductor boundary condition and 
the symmetric solution with vacuum condition on the one hand, and 
the symmetric solution with perfect conductor condition and the 
antisymmetric solution with vacuum condition on the other hand. 

In the following we consider both symmetric solutions using 
the vacuum boundary conditions, as well as antisymmetric ones us- 
ing the perfect conductor boundary condition, which corresponds 
in each case to the most easily excited mode. In the cases where 
we use a vacuum condition we shall sometimes also apply a wind. 
This makes the dynamo somewhat harder to excite and raises C"'* 
from 5.12 to 5.60 for Cu = 0.6, but the associated magnetic he- 
licity flux alleviates catastrophic quenching in the nonlinear case. 
Alternatively, we consider an explicit removal of magnetic helicity 
to alleviate catastrophic quenching. In cases with perfect conductor 
boundary conditions the most easily excited mode is antisymmetric 
about the equator, which corresponds to a boundary condition that 
permits a magnetic helicity flux through the equator. This would 
not be the case for the symmetric solutions. 

3.2 Saturation behavior for different values of Rm 

We now consider the saturated state for a value of Ca that is su- 
percritical for dynamo action. In the following we choose the value 
Ca = 8. The saturation behavior is governed by equation JI7l l. 
Throughout this paper we assume k{/ki = 10 for the scale separa- 
tion ratio. This corresponds to the value 5 in IShukurovet"Z|j2006l) . 
where ki was defined differently. The dynamo saturates by build- 
ing up negative qm when aK is positive. This diminishes the total 
Q in equation JI3t and saturates the dynamo. The strength of this 
quenching can be alleviated by magnetic helicity fluxes that lower 
the negative value of qm. 

We plot in Fig. [3] the dependence of the saturation field 
strength Bsat, defined here as the maximum of at the time 

of saturation. To monitor the degree of quenching we also plot in 
Fig. [3] the Rm dependence of the maximum of the negative value 
of aM at the time when the dynamo has saturated and reached a 
steady state. The maximum value of — oa/ is lowered by about 5% 
from 1.8 to 1.7 in units of rytfci (see Fig.O. Finally, we recall that 
for the dynamos considered here both Bx and By oscillate, but 
their relative phase shift is such that B is non-oscillatory. The nor- 
malized cycle frequency, Cj = uj/rjtki/is also plotted in Fig.|3]as a 
function of Rm- It is somewhat surprising that u does not strongly 
depend on Rm. One may have expected that the cycle frequency 
could scale with the inverse resistive time rjkl. On the other hand, 
for oscillat ory ail dynamos the cycle freque ncy is known to scale 
with rjtkf iBlackman & Brandenburd j2002f) . although that value 
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Figure 3. Scaling of the extremal value of qji./, the .saturation field strength 
i?sat. and the cycle frequency w with B,m and either Cu = 0.6 (solid 
lines) or Cu = (dashed lines). 



could decrease if ?7t(S) is strongly quenched. However, simula- 
tions only give evidence for mil d quenching jBrandenburg et all 
l2008l:lKapvla & Brandenbur3l200ft1 . 

There is a dramatic difference between the cases with and 
without magnetic helicity fluxes. For Cu = 0.6 the dynamo 
reaches asymptotic behavior for large values of Rm, while for 
Cu ~ the saturation field strength goes to zero and max(— qm) 
reaches quickly an asymptotic value corresponding to a level of 
quenching that makes the dynamo marginally excited. 



3.3 Helicity fluxes tlirougli the equator 

We have seen in Sect. |3.1| that in the perfect conductor case the an- 
tisymmetric solutions are the most easily excited ones. The bound- 
ary conditions for antisymmetric solutions permit magnetic helic- 
ity transfer through the equator. However, this alone does not suf- 
fice to alleviate catastrophic quenching unless a sufficiently strong 
flux is driven through the equator. A possible candidate for driving 
such a flux would be a diffusive flux driven by the Vqa/ term. In 
Fig.|4]we plot the Rm dependence of max(— qm), -Bsat, and lj for 
Ka ~ 0.05 and 0. Again, catastrophic a quenching is alleviated 
by the action of a magnetic helicity flux, but this time it is through 
the equator. The maximum value of — qm is lowered by 15% from 
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Figure 4. Same as Fig. [5] but for antisymmetric solutions in a model with 
perfect conductor boundary conditions with Cu = and 
0.05 (solid lines) or (dashed lines). 



2.35 to 2. 15 in units of rjtk\ (see Fig.|4ll. Again, the cycle frequency 
is not changed significantly. 

In Fig.|5]we compare the profiles of ftm, hf, Fm, and Ff for 
the most easily excited solution with vacuum and perfect conductor 
boundary conditions on 2: = L. In all cases we have hm = h{ — 
at the midplane due to symmetry reasons, and at z — L we have 
hm = and hf 7^ 0. It turns out that the magnetic helicity flux of 
the small-scale field is balanced nearly exactly by that of the mean 
field. T his agrees with the expectation of lSlackman & Brandenburg 
(l2003h who argued that both should be shed at nearly the same rate. 

The ad hoc assumption of a turbulent magnetic helicity flux is 
plaus ible and has of course been made in the past jKleeorin et al.l 
but its effect in alleviating catastrophic quenching has not 

S et been seen in earlier three-dimensional tu rbulen ce simulations 
Brandenburg & Dobleijr200ll : lBrandenburgH2001bl) . However, ex- 
cept for the effects of boundaries, the conditions in those simula- 
tion were essentially homogeneous and the gradients of magnetic 
helicity density may have been just too small. It would therefore be 
important to reconsider the question of diffusive helicity fluxes in 
future simulations of inhomogeneous helical turbulence. 
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Figure 5. Mean magnetic helicity densities of mean and fluctuating fields 
as well as mean magnetic helicity fluxes of mean and fluctuating fields as 
functions of z for the S solution with vacuum boundary condition and ad- 
vective flux with Cjj = 0.6 (upper two panels) and for the A solution with 
perfect conductor boundary condition and diffusive flux with ka = 0.05 
(lower two panels). The profiles of hi have been scaled by a factor of 
10 to make them more clearly visible. In all cases we used Ca = 8 and 
Rm = 105. 



3.4 Occasional removal of ht 

Catastrophic quenching can also be alleviated by the artificial 
removal of small-scale magnetic fields. This was first demon- 
strated using three-dimensional simulations of helical turbulence 
( [Brandenburg. Dobler. & Subramanianl 1200 2). In that study the 
small-scale magnetic field was isolated in Fourier space and re- 
moved in regular time intervals. Here we do the same by just reset- 
ting hi to a reduced value, as described in equation il9\ . 

We consider the saturation strength of the magnetic field, 
_Bsat, to characterize the alleviating effect of small-scale magnetic 
helicity losses. Not surprisingly, the dynamo becomes stronger 
(Bsat increases) when the fraction of small-scale field removal e 
is increased (upper panel of Fig. [SJ or the time interval of field 
removal is decreased (lower panel of Fig. |6]l. These dependencies 
follow approximate power laws, 
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Figure 6. Saturation field strengtli versus e for Air^tfc^ = 0.25 (upper 
panel) and versus Atijtk^ for e = 0.1 (lower panel) in a model with Ca = 
8, _Rm = 10^, and Cu = Cs = Ka = 0. 

Bsat/Beq « 0.17e'/^ ^ 0.024 (Afj^t fc? )" (26) 

suggesting that even relatively small amounts of magnetic helicity 
removal in long intervals can have an effect. 

We have also performed some numerical experiments where 
the magnetic helicity associated with the small-scale field is only 
removed near the surface layers. However, in those cases the catas- 
trophic quenching was not noticeably alleviated. This can be ex- 
plained by noticing that, in the absence of additional magnetic he- 
licity fluxes in the interior, there is still a build-up of h{ in the inte- 
rior which quenches the a effect catastrophically. 
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Figure 7. Space-time diagrams for Bx and By, as well as the magnetic 
helicity densities hm and h[ for L = AH, Ca = 8, Cjj = 0.6, and 
Hif oo. The white horizontal line marks the location z = H. Light 
(yellow) shades indicate positive values and dark (blue) shades indicate neg- 
ative values. 



3.5 Magnetic helicity density and flux profiles 

In an attempt to understand further the evolution of magnetic helic- 
ity we have performed calculations where the magnetic helicity flux 
of the fluctuating field was forced to vanish near the surface. This 
was done by choosing a profile for U that goes to zero. However, 
this invariably led to numerical problems. In order to clarify the ori- 
gin of these numerical problems we chose to adopt a taller domain 
L = AH using the profiles shown in Fig.[T] The result is shown in 
Fig. |7] In that case the flux is still able to carry magnetic helicity 
away from the dynamo-active region into the outer layers z > H. 
The cyclic dynamo < z < H operates otherwise very much 
like before in the kinematic regime (Fig.|2]l. Dynamo action is pos- 
sible for Ca > C"'' = 4.32. However, a problem still arises when 
a parcel of positive magnetic helicity that is shed early on from the 
dynamo-active region reaches the upper layers at z « 3H, through 
which now no magnetic helicity can be transmitted. Positive mag- 
netic helicity builds up near z ~ 3-ff until wiggles appear (Fig. (8)1. 
This demonstrates that, once a magnetic helicity flux is initiated, 
there is no way to stop it locally. Instead, wiggles develop, so the 
solution becomes numerically invalid at that moment. 

The fact that positive magnetic helicity is produced is some- 
what unexpected, because for a > the magnetic helicity produc- 



tion is positive definite. However, this can be traced back to the 
term rjtJ ■ B, which is part of f • i? on the right hand side of equa- 
tion Since J ■ B is positive for positive ok, it is clear that this 
term produces positive h{ just outside the range where ok is finite 
and where it would produce h{ of opposite sign. 

In another experiment we adopt a profile for U such that Hu 
is changed from oo to 3H only after a time trjtkl = 25, which 
is when the positive bump of hf has left the domain. The result is 
shown in Fig.|9] Now it is indeed negative magnetic helicity that the 
dynamo tries to shed. However, even though the flux is relatively 
weak, the blockage at z = 3H leads eventually to a problem and 
leads, again, to wiggles indicating that the solution is numerically 
not valid. 

These results suggest that the magnetic helicity flux must be 
allowed to continue through the rest of the domain. Of course, in 
reality there is the possibility of various fluxes, including diffusive 
fluxes that have not been included in this particular model. We note, 
however, that model calculations with finite in equation ( 12 It 
confirm that then Ff{z, t) becomes constant in the outer parts. 
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Figure 8. Evolution of h{ in a model similar to that of Fig. [T] but for 
Hjj — > oo = 3H, so the flux of magnetic helicity of the fluctuating field 
is blocked at z = 3H. Note the emergence of a shock that can eventually 
no longer be resolved and leads to wiggles. The white horizontal line marks 
the location z = H. 




Figure 9. Similar to Fig. [8] but this time the value of H has been changed 
from infinity to Hjj = 3H only at the time rjtk'^t = 25 (marked by a 
vertical white line). Soon after that time the solution develops wiggles that 
invalidate the results after that point. 




Figure 10. Dependence of FmiH) and F{{H) on the shear parameter 
for the S solution in a model with vacuum boundary condition, Ca = 8, 
Cu = 0.6, and Rm = lO^. 



3.6 Magnetic helicity with shear 

It is remarkable that the magnetic helicity fluxes of the mean and 
fluctuating fields were always equally strong and of opposite sign. 
The point of this section is to underline that this is a particular 
property of the dynamo, and would not carry over to afl dy- 
namos. In Fig. [To] we show the fluxes of the model with Ca ~ 8 
and Cu =0.6 where we have varied Cg in the range from —8 to 
+8. 

It turns out that the pre sence of shear gives rise to an addi- 
tional magnetic helicity flux jBerger & RuzmaikinlllOOOh and that 
the perfect correspondence between magnetic helicity fluxes of op- 
posite sign for mean and fluctuating fields is then broken. This ad- 
ditional flux of magnetic helicity is associated with the mean field 
and does therefore not, on its own, alleviate catastrophic quench- 
ing. However, in this model we have neglected magnetic helicity 
fluxes that would be associated with the fluctuating field. An exam- 
ple is the Vishniac-Cho flux whose effect in a mean-field model was 
alread y studied in an earlier paper ( Brandenburg & Subramaniar] 
l2005l) . For C < —2, the oscillating solutions are no longer pre- 
ferred and a new solution branch emerges, where the solutions are 
now non-oscillatory. T hose are also the type of solutions studied by 
Shuk urov et al.l j2006h . where Cs ~ —8 was chosen, correspond- 
ing to the value —2 in their normalization. 



4 CONCLUSIONS 

The present simulations have confirmed that in finite domains mag- 
netic helicity losses through local fluxes are able to alleviate catas- 
trophic quenching. Without such fluxes the energy of the mean field 
goes to zero in the limit of large Rm, while in the presence of 
such fluxes \B\ reaches values that are about 5% of the equipar- 
tition value. We emphasize at this point that this applies to the case 
of an dynamo. For an afl dynamo the mean field can reach 
larger values, depending on the amount of shear. For example for 
the model shown in Fig.[TO]the field strength in units of the equipar- 
tition value rises from 5% without shear to about 36% with negative 
shear (Cs ~ —8), while for positive shear it stays around 5%. We 
also emphasize that the difference between the two cases with and 
without helicity fluxes is rather weak for Rm < 10'^, so one really 
has to reach values around Rm < 10^ or Rm < 10'' to make a 
clear difference compared with catastrophic quenching. Such high 
values of Rm are not currently feasible with three-dimensional tur- 
bulence simulations. 

The other surprising result is that it is not possible to dissipate 
magnetic helicity flux locally once it is initiated. If the magnetic 
helicity flux of the small-scale field has already left the dynamo- 
active domain, it has to stay constant in the steady state. By adding 
a diffusive flux the shock-like structure in the magnetic helicity of 
the small-scale field could be smoothed out, but this contribution 
would then carry the same amount of energy as before, although 
now by other means. 

In the presence of shear there are additional contributions 
to the magnetic helicity flux associated with the mean magnetic 
field. There are first of all the fluxes associated with the mean 
field itself, but those fluxes cannot contribute to alleviating catas- 
trophic quenching on their own. However, earlier work has shown 
that in the presence of shear there are al so additional contribu- 
tions associated with the fluctuating field jVishniac & Choll200ll; 
ISubramanian & Brandenburdl2004 l2006h . Those terms have not 
been included in the prese nt work, because they have already b een 
studied in an earlier paper jBrandenburg & Subramani^l2005l) . 
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Several new issues have emerged from the present study. 
The fact that diffusive magnetic helicity fluxes through the equa- 
tor can alleviate catastrophic quenching is not surprising as such, 
but its effects in alleviating catastrophic saturation behavior in 
three-dimensional turbulence simulations has not ye t been reported 
terandenburg & DobleillioOll : lBrandenburg||2001bl) . On the other 
hand, simulations of forced turbulence in spherical shells with 
an equator did sh ow near-equipartition strength saturations fields 
jMitra et al.ll2009h . although the values of were typically be- 
low 20, so it was not possible to draw conclusions about catas- 
trophic quenching. A new dedicated attempt in that direction would 
be worthwhile using again driven turbulence, but now with a linear 
gradient of its intensity and in cartesian geometry. 

In view of applications to the Sun and other stars, another im- 
portant development would be to extend the present work to spher- 
ical domains. Agai n, some work i n that direction was already re- 
ported in , Brandenburg et al.l | |2007|) . but none of these models used 
diffusive fluxes, nor has any attempt been made to model the Sun. 
This would now be an important target for future research. 
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